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SI.  INTRODUCTION 

The  present  paper  is  concerned  with  propagation  of  chaos  problems  for 
systems  with  an  infinite  number  of  degrees  of  freedom  such  as  strings  or 
spatially  extended  neurons.  The  investigation  of  the  asymptotic  behavior  of 
the  voltage  (membrane)  potentials  of  large  assemblages  of  interacting  neurons 
leads  to  precisely  such  problems  and  provided  the  immediate  motivation  for  the 
work.  Another  example  to  which  the  approach  of  the  present  paper  could  be 
applied  (we  believe)  is  the  Ginsburg-Landau  model  in  hydrodynamics  recently 
studied  by  T.  Funaki  [4]. 

Sections  2  and  3  are  of  an  introductory  nature.  Basic  properties  of  duals 
of  nuclear  spaces  (denoted  throughout  by  <f> '  ,  the  strong  dual  of  a  countably 
Hilbertian  nuclear  space  $)  are  briefly  discussed  and  the  results  of  Kallianpur 
et  al.  [8]  on  the  existence  and  uniqueness  of  the  solution  to  (the  martingale 
problem  posed  by)  a  ^’-valued  stochastic  differential  equation  (SDE)  is 
extended  to  a  system  of  such  equations.  The  principal  results  in  which  the 
infinite  dimensionality  of  our  problem  call  for  special  arguments  are  derived 
in  Sections  3.  4,  and  5. 

In  Theorem  4.1,  the  weak  compactness  of  the  sequence  of  empirical  measures 
1  n 

p  (o> . • )  -—26  (■)  is  established  and  it  is  shown  in  Section  6  (Theorem 

n  n  J=1  X*J(..U) 

6.1)  that  Tjn,  the  law  of  pn(w,»)  converges  weakly  to  the  unique  solution  of  the 
McKean- Vlasov  equation. 

The  infinite  dimensional  (nucelar  space-valued)  version  of  the 
McKean-Vlasov  SDE  is  introduced  in  Section  5.  The  existence  and  uniqueness  of 
solution  of  this  equation  is  investigated  in  detail  in  Baldwin  et  al.  [1].  In 
view  of  the  importance  of  this  result  for  the  propagation  of  chaos,  a  slightly 
different  proof  (with  a  somewhat  stronger  conclusion)  is  given  for  the  special 
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choice  of  the  interaction  term  for  our  problem.  The  main  results  on  the 
propagation  of  chaos  are  given  in  Theorems  6.2  and  6.3.  The  existence  of  a 
unique  solution  to  the  martingale  problem  posed  by  the  system  (3.2.1)  and 
Theorem  5.1  on  the  McKean-Valsov  equation  are  the  key  steps  that  enable 
Sznitman’s  technique  for  finite  dimensional  SDE's  to  be  used  for  the  nuclear 
space  valued  case. 

The  application,  alluded  to  above,  to  the  voltage  potentials  of 
interacting,  spatially  extended  neurons  is  considered  in  Section  7.  For 
reasons  of  space  we  have  limited  ourselves  to  the  mathematics  of  the  problem 
and  excluded  any  discussion  of  the  neurophysiological  implications. 

In  Section  8  we  introduce  the  assumption  that  the  initial  measure  of  the 
system  (3.2.1)  is  pQ-chaotic  and  show  that  the  results  of  the  previous  sections 
hold  under  this  more  general  condition.  This  is  of  importance  in  application 
since  it  is  more  reasonable  to  assume  (as  in  the  case  of  the  neurons)  that  the 

random  variables  Xa(t) . Xa(t)  are  exchan&eable  that  they  are 

identically  distributed. 

It  is  worth  remarking  that  our  results  contain  the  finite  dimensional 
results  as  a  particular  case  and  their  relationship  with  other  available 
results  (e,g.  Sznitman  [16])  is  also  briefly  commented  upon. 

An  cutstanding  problem,  to  which  we  hope  to  return  in  a  later  paper,  is 
that  of  proving  a  fluctuation  or  central  limit  theorem.  The  difficulties  that 
lie  ahead  are  foreshadowed  in  a  recent  paper  by  Kallianpur  and  Mitoma  [7]  that 
establishes  such  a  result  under  restrictive  conditions. 

§2.  PRELIMINARIES  ON  NUCLEAR  SPACES  AND  <t> ' -VALUED  SDE’s. 

In  this  section  we  provide  the  basics  on  nuclear  spaces  and  on  stochastic 
processes  and  integrals  taking  values  in  duals  of  nuclear  spaces  followed  by 
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the  results  of  Kallianpur,  Mitoma  and  Wolpert  [8]  on  the  existence  and 
uniqueness  of  solutions  of  SDE’s. 

2.1  Nuclear  spaces 

Let  $  be  a  real  linear  space  whose  topology  is  given  by  an  increasing 

sequence  11*11  .  r=l,2,...  of  Hilbertian  norms.  Let  4>r  be  the  completion  of  $ 

with  respect  to  11*11  .  Then  is  called  a  countably  Hilbertian  nuclear  space 

(CHNS)  if  the  following  two  conditions  are  satisfied: 

00 

(i)  <t>  =  n  $ 

,  r 
r=l 

(ii)  For  each  r,  there  exists  an  m>r  such  that  the  canonical  embedding  <P  C 

v  '  m  r 

is  Hilbert-Schmidt. 

Let  <f>'  denote  the  strong  dual  of  $  whose  topology  is  given  by  the 
following  family  of  semi-norms: 

|f  L  =  sup  |f(x)  j  where  B  C  4>  is  a  bounded  set  in  <f>. 

x€B 

00 

It  is  well  known  that  $'  =  U  <P  where  $  is  the  dual  of  <P  .  Besides,  the 

.  -r  -r  r 

r=l 

strong  topology  on  <P'  coincides  with  the  inductive  limit  topology  induced  by 
the  canonical  embeddings  <f>_r  C  <P'  .  Let  11*11  denote  the  norm  in  <3>_r .  If  j'r 
denotes  the  canonical  mapping  of  #  onto  its  dual  ^_r-  then  for  u  €  <P_r  and 
*  €  *r. 

uO]  =  <u.Jr«>_r  =  <j_r.u.<A>r 

where  <  •  >  denotes  the  inner  product  in  the  appropriate  space. 

T 

For  any  T>0,  C^,  denotes  the  space  of  all  continuous  functions  from  [O.T] 

to  <P' .  If  {|*la:  a  €  A}  is  the  set  of  semi-norms  defining  the  strong  topology 

T  T 

of  then  by  defining  mxw  =  sup  |x  |  ,  x€C*.  ,  the  space  C.,  is  seen  as  a 

“  OitiT  1  a  9 
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completely  regular  topological  space  under  the  projective  limit  topology  of 

a  £  A} .  ,  denotes  the  space  of  all  continuous  functions  from  [0,°°)  to 

T 

<t>'  .  C.  is  the  Banach  space  with  the  uniform  topolc3y,  consisting  of  all 
-i 

continuous  functions  from  [0,T]  to  4>_£. 

2.2  <ft' -valued  processes. 

Let  (Q,?, (?t) t^Q. P)  be  a  stochastic  basis  satisfying  the  usual  hypotheses. 

Definition "•  An  adapted  -valued  stochastic  process{Mt) is  called  a 

martingale  with  respect  to  (?t)  if  for  each  $  €  <P,  {Mt[$]}^Q  is  a  real-valued 

o  2 

(?t)  martingale.  {M^}  is  called  an  L  -martingale  if  <«>  for  all  t>0  and 

For  a  detailed  discussion  of  ^'-valued  martingales  and  their  properties  we 
refer  the  reader  to  [6]  and  [12]. 

Def inition:  A  continuous  -valued  process  {^t}t>o  is  a  Wiener  process 

with  covariance  Q,  if  the  following  conditions  are  satisfied: 

(i)  WQ=0  a. s . 

(ii)  {Wt[$]}  is  a  one-dimensional  Wiener  process  with  variance  parameter 
Q(*.<fr) .  where  Q(*,*)  is  a  continuous  positive  definite  symmetric  bilinear  form 
on  <£. 

A  result  of  Mitoma  [13]  implies  that  any  -valued  Wiener  process  W  has 
paths  that  lie  in  the  Banach  space  for  some  q<°°,  and  which  are  continuous 

-q 

in  the  $_^-topology  P-a.s.  The  choice  of  q  depends  only  on  the  covariance  form 
Q.  Let  r  2  q  be  a  fixed  integer.  An  important  property  of  the  quadratic  form 
Q  is  that  it  admits  a  unique  continuous  extension  to  a  nuclear  form  on  #r  and 

Q[*.*]  =  (*.Qr*)r  =  (Q^.Q^)r 


(2.2.1) 


5 


for  a  unique  non-negative  trace-class  operator  Qr  on  4*^. 

* 

on  (or,  equivalently,  of  Qr)  is  given  by 


The  trace  norm  of  Q 


lQl.r..r=  achj.hj] 


(2.2.2) 


where  {h^}  is  a  CONS  for  4>  . 

\  j/  r 


2.3  Stochastic  integrals  in  <ft' 


Let  {W^.}  be  a  ^'-valued  Wiener  process  with  covariance  form  Q(*.*)  and  let 

L(4>’,4>')  be  the  space  of  all  continuous  linear  operators  from  4>’  to  4>' . 

2 

For  each  T>0  and  ^  €  <P,  let  denote  the  space  of  progressively 

measurable  processes  H:  IR+xfi  — »L($',$‘)  for  which  EJq  Q[H*^,H*$]ds  <  °°,  where 

H  is  the  operator  dual  to  H  . 
s  s 

Def inition:  The  stochastic  integral  I*?  :  =  dW  (0£s£T)  is  a  <£’-valued 

t  US  S 

2  H  t 

L  -martingale  with  the  quadratic  variation  process  as  <1  Qjj  [♦.4']ds 

s 

where  0,4/]  =  Q[H*$,H*4/]. 
s 

H  T 

There  exists  an  £>0.  depending  on  H  and  T  such  that  I  €  C,  a.s..  If 

v-e 

{hj}  C  $  is  any  CONS  in 


(2.3.1) 


where  the  right  hand  side  is  an  L  -convergent  series  of  I to  integrals. 
Besides, 


t  e , 


>tO.*]  =  2  Sq  (Hs^.hpg(Hs^.hpds  Q[hJ,hp 
i ,  j~l 


(2.3.2) 


2.4  <E‘ -valued  SDE*s. 


We  give  below  the  result  of  Kallianpur,  Mitoma  and  Wolpert  [8]  on  the 
existence  and  uniqueness  of  solutions  of  stochastic  differential  equations. 

For  a  probability  measure  p q  on  4>'  and  a  pair  of  functions 
A:  R+x<f>'  — ►  <t>'  and  B:  K+x<I)'  — ►  L(<f>'  :<£’  )  ,  consider  the  following  SDE: 


dXt  =  A(t,Xt)dt  +  B(t.X  )dW  ] 

xo  =  x<°>  j 


(2.4.1) 


where  Xq  is  a  <f>'-valued  random  variable  with  law  of  X^  given  by  p^,  and  W  is  a 
^'-valued  Wiener  process  with  covariance  form  Q. 

Let  consist  of  all  functions  f:  — ►  P  with  f(u)  =  f(u[$])  where 

f  €  C^(R)  and  $  €  <J>.  The  operator  Lg  is  defined  as  follows:  For  each 

f  € 

Lsf(u)  =  f’(u[>])A(s.U)0]+  ^’(uWlQfB^s.uJfB^s.u)*) 

where  B  (s,u)  :  <P  — »  4>  is  the  adjoint  of  B(s,u),  i.e.,  for  all  v  €  <J>'  and 
<t>  €  <J>. 

v[B*(s,u)$]  =  B(s,u)v[$]. 

Let  fl  =  C. ,  ,  3\  =  Borel  o-algebra  of  cl,,  and  J  =  V  f  and  let 

t>0 

:=  C^, — ►  <t>'  >e  the  canonical  process  defined  by  ir^y  =  y(t)  for  all  y€C^,  . 

If  p€w(cj  ),  then  pTr(.1(A)=p(y  €  cj  :  yf  A),  A  €  3(<f>_k). 

-k  -k 

Definition:  A  solution  to  the  martingale  problem  posed  by  (2.4.1)  is  a 

2 

probability  measure  p  on  C^.  such  that  for  any  f  €  2)^(<l>'),  the  real-valued 
process  =  f(x£)  -  f(Xg)  “  Jq  Lgf(xg)ds  is  a  (n,?,(?t),p)  martingale  with 

=  »0. 

The  following  conditions  are  imposed  on  the  space  <f>,  measure  p^  and  the 


coefficients  A  and  B. 
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Let  (h.)  be  a  CONS  in  <t>  obtained  by  the  Gram- Schmidt  process  applied  to  a 
J  m 

countable  dense  set  (f  .}  in  $.  For  every  j,  we  then  have 

<1 


n . 


f  4  *  "j 


(2.4.2) 


where  n  (depending  on  m  and  j)  <  j  and  llp.ll^  =  0-  Our  assumption  is  the 
J  1 

following: 

(A)  For  each  m  and  p,  (p>m) ,  in  (2.4.2)  lln.ll  =  0. 

J  P 

Let  T>0  be  fixed.  Then,  for  each  sufficiently  large  m>r  where  r  is 
introduced  in  (2.2.1),  there  exists  a  number  0>O  and  an  index  p>m  such  that  for 
all  s , t  i  T, 


(IC)  Initial  Condition:  r  :  =  J  ( 1+ Jlu m) [ ^n ( 3+llullf m) ]2^q ( du) <«» 

(GC)  Coercivity  Condition:  For  each  u  €  j  <t>, 

m 

2A(t.u)[j  u]  +  |Q„,,  J  £  0(l+llu!l2  ) 

v  -m  J  1  B(t,u) v  -nr 

where  j  denotes  the  canonical  map  from  to  <t>  .with  j  as  its 
m  m  -m  — m 


m 

inverse. 


-m 


(LG)  Linear  Growth  Condition-'  If  u  €  $  ,  then  A(t,u)  €  $  and 

v  ’  -m  -p 

IIA(t.u)ll2p  ^  0(l+llull2m) 

|CL,,  a  I  1  0(  1+llull2  ) 

'  B(t,u)  -m,-m  v  -nr 

(JC)  Joint  Continuity  Condition:  A  and  B  are  each  jointly  continuous. 
Further , 

(i)  B(s,u)(v)  €  <t>  if  u.v  €  <P  and 

(ii)  CL.,  -,(♦.♦)  is  continuous  in  u  on  <t>'  for  each  $  €  <t>. 
v  J  B(s,u)v  ’ 


The  following  condition  will  be  needed  in  the  proof  of  uniqueness. 
(MC)  Monotonicity  Condition:  For  all  u.v  €  (C  4>_p) 

(A(t.u)  -  A(t.v).  u-v)_p  +  |QB(t>u)_B(t  v)  Lp.-p  *  0,lu-v,,?P 
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We  give  below  the  main  result  of  Kallianpur  et  al.  [8], 

Theorem  2.4.1.  Assuming  conditions  (A),  (IC),  (OC) ,  (LG)  and  (JC),  there  exists 
a  weak  solution  to  the  stochastic  differential  (2.4.1).  Besides,  it  has  the 
pathwise  uniqueness  property  if  ((MC)  is  satisfied. 

Next,  we  give  a  moment  bound,  followed  by  a  tightness  result  both  of  which 
are  due  to  Baldwin  et  al  [1]. 

2k 

Theorem  2.4.2.  Let  k£l  and  EIIXqII^  £  c^  <  «.  Then,  under  all  the  conditions 
of  Theorem  2.4.1, 

E  sup  NX  11^  £  (2C  +  l)exp((136k2  -  4k)0T)-l 
O^T  s  m  k 

Remark  2,4,1:  EIIX^I2^  *  (2(^+1  )e4k^k~1  ^0t  -  1  for  each  0<t£T. 

Theorem  2.4.3:  Assume  that  the  coefficients  associated  with  the  equations 

X  t  =  X0  +  Sofis.tyds  +  J*Bn(s,x")dWj 

and 

Xt  =  X0  +  ^0A(S,Xs)ds  +  JoB(s.Xs)dWs 
satisfy  the  conditions'. 

1)  Conditions  (IC),  (CC) ,  (LG),  (JC)  and  (MC)  hold  as  stated  where  the 
constants  and  indices  are  independent  of  n. 

2)  Xn  X„ 

3)  If  and  Q  denote  the  covariance  forms  of  {W^}  and  (Wg)  respectively,  then 
Qn  converges  to  Q. 

4)  For  each  s£[0,T]  and  <p  €  <P,  An(s,*)[£]  converges  continuously  to  A(s,  •}[$]. 

5)  For  each  s€[0.T]  and  ^  €  <J>,  (Bn(s,*))*£  converges  continuously  to  B*(s,  *)<f>. 

Then  P* (Xn)_1  =>  P*X_1  in  cj 

-P 
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Remark  2.4.2:  If  T>0  fixed,  the  solution  of  SDE  (2.4.1)  namely  X  will  have 
T 

paths  in  a.s.  where  p  is  the  index  that  appears  in  the  conditions. 

~P 

Remark  2.4.3.  Throughout  the  paper,  the  notation  5-»  is  used  to  denote 

convergence  in  distribution  of  random  variables  whereas  the  notation  =>  is  used 

to  denote  weak  convergence  of  measures.  Thus,  x  is  equivalent  to 

P*(X  )  1  =>  P*X  1,  where  PX  ^  and  PX  *  denote  the  law  of  X  and  X  respectively, 
n  n  n 

We  adhere  to  this  notation  even  when  the  random  variables  are  measure-valued. 

3.  SYSTEMS  OF  <t> '  -VALUED  SDE  ’  s 

The  aim  of  this  section  is  to  extend  the  results  of  the  previous  section 
to  a  system  of  stochastic  differential  equations  which  is  done  by  first 
introducing  the  Cartesian  product  of  nuclear  spaces. 

3. 1  Cartesian  product  of  nuclear  spaces. 

Let  <t>  denote  the  nuclear  space  introduced  in  Section  2.  Consider  the 
linear  space  <f>x$  with  coordinatewise  linear  operations.  Let 

"(VV?  =  "Vr  +  N*2Nr  for  (3.1.1) 

An  increasing  sequence  of  Hilbertian  norms  is  thus  defined  on  $x<$  which 

preserves  nuclearity  of  $x$.  To  see  this,  let  $rxr  be  the  completion  of  <£x<f>  in 

the  product  r-norm  given  by  (3.1.1)  for  all  r>l.  Clearly,  <*’rxr=  #rx$r  and 

=  fl  (♦  x$  ),  Given  n>0,  if  m>n  such  that  the  canonical  injection 
rll 

i:  $  C  <t>  is  Hilbert-Schmidt,  then  the  injection  C  ($x<£)  is  also 

m  n  J  v  /m  v  •'n 

Hilbert-Schmidt. 

Let  (<t>x<t>) '  denote  the  strong  dual  of  4>x4>  so  that  (<J>x<!>)  ’  =  U  (<t>  x<P  )  1  .  If 

rll  r  r 

2  €  ($rx<J>r)'  ,  we  can  uniquely  determine  two  linear  functionals  2^  and  2^  in  $  r 


such  that 
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=  V*l)  +  W  for  a11  *1  •  *2€  *r  (3.1.2) 

Likewise,  given  £.  and  £0  in  i>_  ,  there  exists  a  unique  £  in  (4>  x<£  ) 1  such  that 
x  z  -r  r  r 

(3.1.2)  is  satisfied.  In  short,  there  exists  an  isomorphism  between  ($  *$r)' 

and  #  x<f>  which  is  written  as  (4>  xtf>  )  ‘  =  <I>  x  <f> 

— r  —  r  r  r  —  r  —  r 

The  above  isomorphism  is  not  just  algebraic  but  topological  as  well,  if 

4>_rx  4>_r  is  equipped  with  the  product  -r  norm.  i.e. 

Il(u1xu0)ll^  slluJI^  +  llujl^  .  We  thus  get 
x  z  “i*  x  r*  z  r 

($x<f>)'  =  U  ($  x4>  )’  S  U  ($>_  x  $  )  (3.1.3) 

r*l  r  r  r^l  r  "r 

Besides,  U  ( <t>_  x  <J>  )  =  $'x  set-theoretically.  To  see  the  topological 

r*l  r  "r 

equivalence,  consider  a  neighbourhood  of  zero  in  i.e.  Let  and  be 

two  bounded  sets  in  $,  and  e>0  be  given.  Consider  the  set 


A  =  ((£,,£0):  sup  <  e,  sup  |*9(*)  <  e)  C  $'x 


For  any  fixed  r£l,  3  af  Aj  C  l)^llr  <  ar)  1=1,2  so  that 


A2((VV:  "V-r  "V-r  < 


Thus 

n2((«1,«2)--  ll«1M?  +  <  fc2/ai> 

On  the  other  hand, 


(3.1.4) 


AC  {(««):  sup  \£  (*  )+«  (* 9)|  <  2e}  (3.1.5) 

(*r*2)  €  AixA2 


and  AjX  Ag  is  a  bounded  set  in  $  x  <f>.  (3.1.4)  and  (3.1.5)  give  us  the 
topological  equivalence  of  U  ( <t>_  x<J>  )  and  tf>'x$’.  Therefore  (3.1.3)  implies 

r*l  r  "r 


that 
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(4x4)'=  U  (*  x  *  )'  “  U  (4  x  4  )  =  4'x  4*.  (3.1.6) 

r^l  r  r  r^l  ~r  ~r 


Equation  (3.1.6)  carries  over  for  any  finite  number  of  Cartesian  products. 


3.2  System  of  SDE’s 

Let  Xj  be  4 '-valued  processes,  l£j£n  governed  by  the  SDE 

dXj( t)  =  (a(t,Xj(t))  +  I  2  b(t.tfj(t).xj(t)))dt 


+  (®(t.xj(t))  +  i  2  c(t.xj(t).x£(t)))dwj 


(3.2.1) 


and  x”(0)  l^j<n  being  iid  ^'-valued  r.v.'s  with  law  of  Xj(0)  given  by  the 
probability  measure  {W^},  l£j£n.  are  independent  copies  of  a  Brownian 

motion  with  Q  as  the  covariance  form.  Besides, 


a:  IR+  x  4'  — »  4' 
b:  R+  x  4'  x  4'  — *  4' 
a :  K+  x  4'  — »  L(4' :  4’) 
c:  IR+  x  4'  x  4'  — »  L(4' :  4') 


Let  (X^( t) , . . . ,X^( t))  €  4'x...x  4'  be  a  solution  of  (3.2.1).  Then  the 
isomorphism  given  by  (3.1.6)  between  4'x...x  4'  and  (4  x...x  4)'  we  can  write 
the  system  of  SDE’s  (3.2.1)  as  follows: 


dX^  =  (a(t.xj)  +  b(t.x"))dt  +  (a(t,X")  +  c(t.xJ))dWt  (3.2.2) 

where  initial  value  Xq  is  a  (4x...x  4)'-valued  r.v.  such  that  Xq  is  isomorphic 
to  (Xj(0) . XJJ(O))  €  (4'x. . .x4' ) ,  and  W  is  the  (4  x...x  4) '-valued  Wiener 


process  described  below. 
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Wt(r)  =  2W {(f). 

1=1  c  J 


By  the  independence  of 


n 


<W>|>.*]  =  t  2 


j=l 


j'rjJ 


Besides.  {X^}  is  a  (<P  x...x  <$)'-valued  process.  The  coefficients  appearing  in 
(3.2.2)  are  given  as  follows: 
a:  IR+  x  (<P  x. .  .x  <P ‘ )  — >  (<£  x.  . .  x  <J>)  1 

b:  IR+  x  [<t>  x...x  0)  '  — ►  (<t>  x.  .  .x  <f>) ' 

a ■  IR+  x  ($  x...x  4>)  '  — ►  L((tf>  x...x  $)  *  :  (4>  x...x  $)  ‘  ) 

c":  IR+  x  (<P  x.  .  .x  $}  *  — *  L(($  x.  .  .x  ♦)  * :  ($  x. .  .x  $>)  ’  ) 

For  <#>  =  (<p . . .  .  ,<p  )  and  $=($, . ♦  )  €  <p  x. .  .x  $  and  u  =  (u  . u  )  and 

/v  a  **  *  n  /x/  a 

v  =  (v. . v  )  €  ($  x. .  .x  <t>)  '  ,  we  have 

*v  ^ 


n' 


a(t,u)[<p]  =  2  a(t,u  )[>  ] 

~  ^  J  J 


b(t.u)[<p]  =2^-2  b(t,u  ,u  )[<p  ] 
j=l  i=l  J  1  J 


»v  <\f 


n 


^(t.u){v)[<^]  =  2  ^’UjXVjJOj] 


<V  'V  »v 


n  n 


so  that 


c(t.u)(v)[>]  =  -  2  2  c(t,u  ,u  )(v  )J>  ] 

~  ~  ~  j=l  i=l  J  J  J 

^o(s.u)dwsW  =  ^2^/g  CT(S'Uj)dW^[<Pj] 

Sl  c( s.u)dWsW  =  i  Ijl  e(S.uJ.u1)dwJ[».] 
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Under  the  conditions  given  below,  there  exists  a  weak  solution  Xn  of 
(3.2.2).  Such  a  solution  will  also  be  shown  to  be  pathwise  unique  and  thus 
ensures  a  unique  strong  solution  {X^}  €  (<£  x...x  <£)  *  to  the  equation  (3.2.2) 
(see  [5]). 

Analogous  to  (3.1.6),  the  isomorphism  C^x  x$)'  S  C^>, x.  .  . xC^,  is  easily 

T 

established.  Towards  this,  fix  T>0  and  consider  y€C^  x  x  <P)'  so  t^iat  ^or 

each  0£t£T,  y(t)  €  (<£  x...x  <f)  ’ .  By  using  (3.1.6),  y(t)  is  isomorphic  to,  say, 

(y^t) . y  (t))  €  4>'x...x  for  0  i  t  i  T  and  so,  lim  y(t)  =  y(tQ)  is 

n  wt0 

equivalent  to  lim  y^( *0^  ^or  eac*1  l£k£n.  Let  {|*la:  a€A}  be  the  set  of 
t-t0 

semi-norms  defining  the  strong  topology  of  Set  wy,«  =  sup  |y,(t)|  for 

K  a  0£t<T  K  “ 

T  2  n  2 

y,  €  C. , ,  1  i  k  $  n  and  a  €  A.  Define  *y*  =  2  »y,iH  for  each  a  €  A.  By 

iC  9  CL  i  «  K  CL 

k=l 

replacinglMlr  by  »»«a  in  (3.1.4),  the  arguments  used  in  deriving  (3.1.6)  hold 

T 

in  the  present  context  as  well  and  thus  C^x  x^y  is  isomorphic  to 
T  T 

C^.x.-.x  C^,  .  C^,x...xC^,  and  x  x  $)  ’  e£lulPPet*  with  the  projective  limit 

topologies  of  (C^,x...x  C^,:  T€IN}  and  (C^,x...x  C^:  T  €  IN}  respectively,  are 
therefore  isomorphic. 

Let  X^t),  l^j^n  and  t>0,  be  such  that  for  each  </>  =  (<#>. . f  )  € 

<t>  X...X  <J>.  x”(<p)  S  2X*J  [<p  ].  Then  X"(*)  €  C.  .  for  1  £  j  £  n  and  solve  (3.2.1). 

The  conditions  for  the  existence  and  uniqueness  of  solutions  of  SDE  (3.2.2)  are 
as  follows:  Conditions  (A)  and  (IC)  of  Section  2  are  assumed  to  hold.  It  is 
easy  to  note  that  condition  (A)  on  $  implies  that  on  <f>  x...x  $.  Likewise 
condition  (IC)  of  Section  2  implies 

X  (1+llullf  )[«n(3+llullf  )]2  p  (du)  <  « 

ni  /v  m  o 

($  X.  .  .X  $)  ' 
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where  u€($x. .  .x<£) '  and  Hq  =  S£(Xq).  We  will  call  (A)  and  (IC)  as  (SA)  and  (SIC) 

where  S  stands  for  system  of  stochastic  differential  equations.  The 
coefficients  a,b,c  and  a  are  assumed  to  satisfy  the  following  conditions: 

For  any  T>0,  3  Lj,  \  r  such  that  for  each  m  £  Lp,  3  a  number  0  and  an  index  p 
(note  that  0  and  p  depend  on  m)  such  that: 

(SOC)  For  u.v  £  and  0  i  t  £  T, 

2a(t.u)[j_mu]  +  2b(t.u.v)[j_mu]+ 


m 


|Q  .  ,  J  *0(l+llull2  +P(t.u.v)llvll2  ) 

1  a(t.u)+c(t,u,v) '-m.-nr  v  -m  y  -my 


where 


/3(t.u,v)  =  \ 


0  if  b(t.u,v)  =  c(t.u.v)  =  0 
1  otherwise 


(SLG)  Let  u.v  €  and  0  £  t  i  T.  Then  a(t,u)  €  $_p  and  b(t,u,v)  €  $_p. 

Besides. 


Ha(t,u)ll2p  *  0(l+llullfra) 

ilb(t.u.v)ll2p  1  0(l+llullfm  +  Hvll2m) 

|Q  J  £  0(l+llull2  ) 

1  o(t, u)'-m,-m  -my 

|Q  J  S  0(l+llull2  +ilvll2  ) 

1  c( t.u.v) '-m.-m  v  -m  -m' 


(SJC)  a.b.c  and  a  are  jointly  continuous  functions.  Further, 

(i)  For  u.v.w  €  <P_m  , 
o(t.u)(w)  e  <t>_m 
c(t,u,v)(w)  e  $_m 

<“>  Qo(t.  continuous  in  u  on  <P'  and  u  *s 

continuous  in  u  on  for  each  £  €  <t>. 

The  following  condition  is  needed  to  prove  the  uniqueness  of  solutions. 
(SMC)  For  urvru2.v2  €  $_m  (  C  *_p) 
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(a(t.u1)-a(t.v1).Vv1)_p  +  i  6  “VV 


and 


(b( 


t.uru2)-b(t.vrv2).u1-v1)_p  +  lQc{ t .u2)-c( t.Vi .v2)  Lp,_p 


*  e  {iiUl-v1uf  +  Iiu2-v2n?  } 


2 

-P 


It  is  easily  verified  that  the  above  conditions  will  imply  (OC) ,  (LG), 
(JC)  and  (MC)  for  the  coefficients  a,b,c  and  a  for  equations  (3.2.2),  so  that 
by  Theorem  2.4.1  we  get  existence  and  uniqueness  of  solutions  for  the  SDE 
(3.2.2).  The  moment  bound  given  in  Theorem  2.4.2  becomes 

E  sup  IIX1J(t)ll2k  <  (2C,  + 1 ) exp ( ( 136k2-4k) 0T) - 1  (3.3.1) 

Oit$T  J  “m  k 

n  2k 

where  =  EIIX^ (0) H_m  <  “  and  is  independent  of  n.  In  fact,  the  following 
bound  can  also  be  derived  by  an  obvious  modification  of  the  proof  of  (3.3.1): 
Under  the  Condition  (IC) 

E  sup  HX"(s)-x"(t)ll2  £  c6*  (3.3.2) 

|  t-s |<6  J  j  -p 

0£t,s£T 

where  c  is  a  constant  independent  of  n.  To  see  this,  note  that  if  T  £  t  >  s  > 
0  with  t-s  £  5,  then 

Xj( t)  =  X*J(s)  +  ^{(a(s).X^(s))  +  i  J^bfs.X^sJ.X^sJJJds 

+  sl  (ct(s,Xj(s))  +  i  ^c(s,x"(s).x"(s))}dwJ. 

We  get  (3.3.2)  by  the  familiar  route  namely,  via.  Doob’s  inequality,  Jensen's 
inequality  wherein  the  condition  (SLG)  is  used  crucially. 

Note:  0<T<°>  will  be  kept  fixed  till  the  last  paragraph  in  Section  6.  Thus  X^ 

J 
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l^j^n.  the  solution  of  3.2.1  will  each  have  paths  in  a.s.  where  p  is  the 


~P 


index  that  appears  in  the  conditions. 


§4.  WEAK  COMPACTNESS  OF  EMPIRICAL  MEASURES. 


Let  X?  denote  the  solution  of  the  SDE  (3.2.2)  so  that 

X?  =  (X*(-) . xn(*))  €  C$  *•••*<£  where  x"(*).  l<j<n  solve  the  SDE’s 

n  -p  -p  J 

(3.2.1). 

T  T  T 

Let  28(C^  )  be  the  Borel  a-algebra  of  .  For  u  €  fi,  B€26(C^  ),  define 

"P  ~P  ~P 

the  empirical  measure 


V  (cj.B)  =  —  2  5 
n  n  .  .  vn,  > 

j=l  Xj(».u) 


(B). 


(4.1) 


For  any  k  2  1.  let  ir( C.  )  be  the  space  of  all  probability  measures  on 

-k 

T 

C,  equipped  with  the  topology  of  weak  convergence  of  measures.  Likewise 
-k 
T 

Tr(C^.)  will  be  the  space  of  all  probability  measures  equipped  with  the  topology 

of  weak  convergence  of  measures.  Note  that  the  canonical  injection 

i:7r(cT  )  C  tt( cT  )  is  continuous  if  k  <  B.  To  see  this,  let  X  €  tt(cT  )  such 

-k  -e  n  -k 

that  X  =>  X  in  ir(cT  ).  Therefore,  for  all  f  €  0,(0*  )•  X  f(y)X  (dy)  — » 
n  -k  b  -k  CT  n 

*-k 

T 

X  f(y)X(dy)  as  n — »  If  g  €  C^(C^  ),  let  j:  tf>  4>_^  be  the  continuous 


-e 


-k 


canonical  injection  of  into  so  that  the  composition  g*j  €  ) 

-k 

Besides.  X  (g*j)(y)X  (dy)  =  X  g(y)X  (dy)  and  X  (g* j) (y)X(dy)  = 

T  n  T  n  T 

-k  -e  -k 

T 

X  g(y)X(dy),  so  that  X  =>  X  in  tt(C.  )  as  n  — ►  «. 

„t  n  -e 


-e 
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T 

p  is  a  random  measure  with  Pn(u>*)  €  ir(C^  )  for  each  gj  €  Q  and  n  >  1. 

n  n  -P 

T 

Let  u  (B)  =  Ep  (B)  for  all  B  €  28(C^  ).  Let  r j  •'=  i£(pn) ,  the  law  of  p .  ,  i.e., 
n  n  ~P  n  " 

T 

the  probability  measure  on  irfC^  )  induced  by  the  random  measure  pn-  Thus, 

-P 

r7n€w(w(cj  )) 

~P 

Theorem  4.1:  Under  the  conditions  (SIC),  (SA) ,  (SCC),  (SLG),  (SJC)  and  (SMC), 
we  have 

T 

(a)  the  sequence  {v  }  is  tight  in  v(C^  )  for  some  q  £  p. 

n  -q 

T 

(b)  the  sequence  (q  }  is  weakly  compact  in  ir(ir(C^  )). 

n  -q 


Proof :  (a)  For  any  set  B  €  28(C^,), 


n 


yn(Bw)  =  E*in(B")  =  ^  2  P(Xj  €  B^)  =  P(X*€BC)  since  X“,  l<j<n.  are  identically 

T 

distributed  random  variables.  Therefore,  the  tightness  of  (i>n)  in  ir(C^,)  is 
equivalent  to  the  tightness  of  {X”)n^ ,  i.e.,  of  the  probability  measures  {Pn} 


on  v(C^.)  where  Pn  =  law  of  X^. 

By  a  result  of  Mitoma  [14],  the  tightness  of  {Pn}  is  equivalent  to  the 

tightness  of  {P  it ■  }  on  v(C^)  for  each  9  €  <P,  where  ir^:  C^,  — *  with 

ir  (u)  =  u[*].  The  tightness  of  (P11  tt . ^ }  follows  by  verifying  the  following  two 
9  9 

conditions: 

(i)  Given  a  >  0,  3  a  >  0  such  that 


sup  Pn(y  €  cl,:  sup  |y  [♦]  |  >  a)  <  e 
n  V  OitiT 


(ii)  Given  e  >  0  and  p  >  0,  3  5  >  0  such  that 


sup  Pn(y  €  c£.:  sup  |y  |  -ys[*3 I  >  *)  i  P- 
n  |t-s|<6 
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Condition  (i)  is  verified  by  noting  that  E  sup 

OitiT 


HX?{t)!l?p 


is  finite  and 


independent  of  n.  To  verify  ( i i ) ,  note  that 


p“(y: 


sup  |yt[<»ysi>]  I  ±  e)  i  Pn(y;  sup  lly  -y  II 11*11  ^  £) 

|t-s|<6  |t-s|<6  v  v 


H\C  n  n  9 

— =2-e  sup  nx?(t)  -  x”(S)ir  i  P 

e  |t-s |<6  1  1  P 


9  2 

if  6  is  such  that  E  sup  lix!*(t)  -  x”(s)ir  *  . 

|  t-s  I  <6  1  1  "P  ll*H 

p 

By  the  moment  bound  given  in  (3.3.2)  we  have  that 


E  sup  BX"(t)  -  X”(s)ll2  <  cliA 

| t-s | <6  1  1  -P 

so  that  a  5  as  desired  does  indeed  exist  once  we  are  given  e.p  and  *. 

T 

The  tightness  of  (un)  is  thus  established  on  the  space  Tr(C^,).  For  any 
*  €  *,  and  any  given  e  >  0. 


P(  sup  |x"(t)[*]|  >  e)  S 
0<,t$T  J 

S 

S 

9  2 

if  11*11  i  - E- - Toogt-  Thus  by  Mitoma  ([14];  Remark  (R.l)),  {v  }  are 

m  (2C1+l)e1^01  n 

uniformly  m-continuous  and  hence,  are  uniformly  p-continuous  as  well. 

T 

Therefore,  there  exists  an  index  q  £  p  such  that  {v  }  is  tight  in  tt(C^  ). 

r‘  -q 

(b)  For  the  second  part  of  the  theorem,  let  us  look  upon  Pn(u.*)  as 


P(  sup  HX"(t)ll  11*11  >  e) 


n*11:  „  9 

— E  sup  lix"(t)ir 

s2  onh  J 


— 2~  l2Ci+1)e  S  P 

e 
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T  T  T 

ir(C^  )-valued  random  variables.  Then  €  ir(C^  )  and  €  tt(it(C^  ))  for  all 

-q  n  -q  n  -q 

n  1  1. 

Note  that  un(B)  =  /  X(B)T?n(dX)  V  B  €  3  tt(cJ  ).  (4.2) 

T(Ct  )  n  -0 

-q 


Using  part  (a)  of  the  theorem,  for  each  j  £  1,  there  exists  a  compact  set  K.  in 

T  _  ^ 

C#  such  that  i>n(K*p  =  J*  X(Kj)dqn(X)  <  e/j  where  e  >  0  is  given.  Let 

K={X'-  MKj)  M  -  jV  j)  where  Kj  ’s  can.  WLOG,  be  taken  to  be  increasing  sets. 

T 

K  C  Tr(Cp  )  is  compact  since  closed  tight  subsets  of  probability  measures  on  a 

-q 

T 

complete  separable  metric  space  (in  our  case,  on  )  are  compact  (see  Chapter 

-q 

II,  Theorem  6.7  in  [15]). 


qn(KC)  =  P(pn  €  KC)  1  2  T7n(X:X(K.)  <  l  ~  j) 

j  •“  1 

00 

=  2  n(X:X(K  )  l  I) 

J=1  n  J  J 


f  MKj) 

1  um 


d\ 


i  e  2  ~  <  2e. 

j=l  J 


T 

Thus  tightness  of  {q^}  in  7r(Tr( C ^  ))  ensues. 

-q 

T 

Note  that  ^(C^  )  equipped  with  the  topology  of  weak  convergence  is  a 

-q 

complete  separable  metric  space  (see  Ch.  II,  Theorems  6.2  and  6.5  in  [15]). 
Tight  subsets  of  probability  measures  on  a  complete  separable  metric  space  are 
relatively  compact.  The  proof  is  thus  complete. 


20 


Remark  4.1.  Let  (q  }  be  the  subsequence  given  by  the  above  theorem  so  that 

T 

q  =>  q  (say).  Let  p  denote  the  ir(C.  )-valued  random  variable  whose  law  is 

"k  % 

given  by  rj.  Thus  q  =>  17  is  equivalent  to  saying  that  p  — *  p. 

\  "k 


Remark  4.2.  By  a  well-known  theorem  of  Skorohod,  there  exists  a  probability 

T 

space  on  which  are  defined  irfC^  )-valued  random  variables,  say  {Z  }  and  Z 

-q  He 

with  law  of  Z  =  q  and  law  of  Z  =  tj  such  that  Z  — *  Z  a.s.  .  Using  this 

"k  "k  "k 

representation  and  applying  Fatou's  lemma,  we  get 


E  sup  J*  Hy  11^  Z(dy)  =  X  {  sup  X  lly  II2  X(dy)}q(dX) 

O^slT  pT  S  ~q  ,rT  .  O^s^T  rT  q 

-q  -q  -q 

i  lira  E  sup  J  lly  II2  Z  (dy) 
k-»»  0£s£T  pT  q  "k 

-q 

=  lilS  X  {  sup  X  lly  II2  X(dy)}q  (dX) 

S"  ^ 

-p  -p 

1  ^  "k  2 

$  lim  E  sup  —  2  MX.  (s)ir 

k-*»  0£s<TT  “k  j.l  J 


$  (2C1+l)exp(1326T)  (4.3) 


by  using  (3.3.1).  The  inequality  (4.3)  holds  with  and  0  remaining  the  same 
if  T  is  replaced  by  any  t  on  both  sides  of  the  inequality  (4.3),  0£t£T. 


Remark  4.3:  In  case  p  is  a  degenerate  random  variable,  and  Y  is  a  -valued 

-q 

random  variable  with  i£( Y)  =  p,  then  Remark  (4.2)  implies  that  for  each  0<t<T. 


E  sup  IIY  II2  $  (2C  +l)exp(1320t). 
Oisit  S  q 


(4.4) 
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§5.  THE  McKEAN- VLASOV  EQUATION 

Let  {Y  :  0^t<T)  be  a  ^'-valued  stochastic  process  that  solves  the 
following  SDE  known  as  the  Global  McKean-Vlasov  equation:  For  0  <  t  <  T, 

Yt  =  Y0  +  Sq  A(s.Ys,^(Y))ds  +  B(s.Ys,tf(Y)}dWs  (5.1) 

where 

A(s,u,X) :  [O.T]  x  x  t r(cj .  )  — »  <f>' 

B(s.u.X):  [O.T]  x  <P‘  x  tt(cJ,)  — >  L($'  :  <P') 

and  W  is  a  ^'-valued  Wiener  process.  S£(Y)  denotes  the  law  of  Y,  and  Y^  is  a 
<t>‘  -valued  random  variable. 

The  local  McKean-Vlasov  equation  is  of  the  form: 

Yt  =  Y0  +  ^0  A(s.Ys.i£(Ys))ds  +  Sq  B(s.Ys,S£(Ys))dWs  (5.2) 

for  Olt^T. 

By  uniqueness  of  solutions  of  the  SDE  (5.1).  we  mean  the  following: 

T 

For  each  X  €  ^(C^.),  let.  for  O^t^T, 

^  =  Y0  +  J0  A(s.^.X)ds  +  Sq  B(s.Y^.X)dWs.  (5.3) 

T  Ai  Ao 

Suppose  there  are  X^  and  X^  6  w( C^,)  such  that  X^  =  Sf(Y  )  and  =  if(Y  ). 

Then  X^X^. 

Existence  and  uniqueness  of  solutions  of  equation  (5.1)  in  full  generali  ty 

will  appear  in  Baldwin  et  al.  [1].  Here,  we  content  ourselves  with  the 

T 

following  choice  of  A  and  B:  For  each  X  €  ^(C^,),  let 

A/ 

A(s,u,X)  =  a(s,u)  +  b(s,u,X) 

A/ 

where  b(s.u.X)  =  J  b(s,u,y  )X(dy)  and 
T  S 

c<y 

B(s,u,X)  =  a(s,u)  +  c(s.u,X)  (5.5) 
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where  c(s,u,A)  =  X  c(s,u,y  )X(dy).  Besides,  we  assume  that  there  exists  a 
T  s 

Ci’ 

sufficiently  large  number  M(m) ,  possibly  depending  on  m,  where  m  is  as  in 

subsection  3.3,  such  that  for  each  u.v  €  $  ,  and  0<s£T, 

— m  ”  ~ 


b(s,u,v)  = 


b  (u.v)  if  llu-vll_m  £  M(m) 

A 

bs(u,v)  otherwise 


(5.6) 


where  b  (u.v)  and  b  (u.v)  are  functions  of  u.v  with  tlb  (u,v)ll_  <  C(m)  for  each 
s  s  s  m 

s  €  [O.T].  Likewise 


c(s,u,v)  = 


cs(u.v) 

cs(u.v) 


if  llu-vll  £  M(m) 
-m  ' 

otherwise. 


(5.7) 


with  c  (u.v)  and  c  (u.v)  are  functions  of  u.v  with  |Q^  |_  _  £  C(m)  for 

s  s  /  »  in,  m 

c  (u.v) 
sv  ’ 

each  s  €  [O.T]. 

With  b  and  c  as  above  we  first  note  that  b  and  c  exist  and  are  finite.  To 
see  this,  consider 


b(s,u.X)  =  /  b(s,u,v  )dA(v)  = 


X  b(>s ,u, v)dAir  (v) 
<t>' 


Since  u  €  ,  there  exists  an  index  k  such  that  u  €  and  such  a  k  can  be 

chosen  to  be  sufficiently  large. 

Using  such  an  index  k  in  the  place  of  m  in  our  conditions  given  in 
subsection  (3.3)  as  well  as  (5.6)  and  (5.7),  we  get  for  each  <?  €  <P  that 


|b(s,u,X)[<#>]  |  £  X  |b(s,u,v)[<#>]  |dX7T  (v) 


=  X|b(s,u.v)[<p]  |dXir'  (v) 
A  s 


23 


+  f  |b(s,u,v)[«f]  |dXirs^(v) 
AC 


where  A  =  {v:  llu-vll_k  £  M(k)}.  Continuing, 


i  S  Ml 

A 


Pk 


( l+llull_k+llvll_k)dAir  1  (v)  +  C(k) 


by  using  (SLG)  with  £  k  as  the  index  that  corresponds  to  k.  Continuing, 

<MI  (M(k)  +  1  +  2l!ull_k  +  C(k)  < 

Pk 

Likewise  one  can  establish  the  finiteness  of  c(s,u,X)  by  showing  that 

|Q_  |_k  _k  <  00  whenever  u  €  ^_k- 

c(s,u,X) 

Such  a  choice  of  b  and  c  makes  physical  sense  in  it  that  a  pair  of 
particles  far  apart  interact  boundedly.  This  choice  includes  in  particular  the 
case  where  b(s,u.v)  and  c(s,u,v)  are  both  bounded  in  the  sense  that 


llb(s,u.v)ll_m  1  C 


|Q-  I  m  m  £  c 

,  ■.  -m ,  -m 

C  (u,v) 


for  all  0£s<T,  u  €  <P'  and  v  €  <P' .  We  assume  that  the  functions  a,b,c  and  a 

satisfy  the  conditions  (SCC) ,  (SLG),  (SJC)  and  (SMC).  It  is  then  a  routine 

matter  to  check  that  A  and  B,  as  defined  above,  satisfy  (CC) ,  (LG),  (JC)  and 

(MC)  as  listed  in  Section  2  with  the  same  indices  and  with  constants 

T 

independent  of  the  measures  X  €  irfC^.).  We  need  and  hence  introduce  the 
following  additional  condition: 

(SJC)'  (i)  For  each  X  €  ^(C^.).  b(s,u,X)  and  c(s,u,X)  are  jointly  continuous 

in  s  and  u.  (<P.<P)  is  continuous  in  u  on  <t‘  for  each  <p  €  <f>  and 

c(s,u.X) 

s  €  [0,T] . 
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(ii)  For  each  $  €  $,  u  €  and  s  €  [O.T],  b(s,u,Xn)[$]  converges 

/V 

continuously  to  b(s,u,X)[$]  and  c  (s.u.A^)^  converges  continuously  to 

T 

c  (s,u,X)$  as  =>  X  in  irfC^,}. 


Theorem  5.1:  Consider  the  SDE  (5.3)  with  A  and  B  specified  by  equations  (5.4) 
to  (5.7).  Assume  conditions  (SA),  (SIC),  (SOC),  (SLG) ,  (SJC).  (SMC)  and  (SJC) 1 

4 

and  that  EIIY_ II  where  C  is  a  positive  constant.  Then, 

0  -m 

(i)  the  McKean-Vlasov  equation  defined  by  the  equation  (5.3)  admits  a  solution. 

X  T 

The  solution  Y  lies  in  CA  a.s. 

<P 

-P 


(ii)  The  solution  is  unique  if  the  following  additional  condition  holds 
(MCjt):  For  all  u.v  €  $_m,  and  €  ir(C^,),  0£s£T, 

^(s.u.^)  -  b(s.v.C2),u-v>  +  IIQ  II 

cfs.u.CjJ-cfs.v.^) 


*  Gr{Ks(f1.C2)llu-vll  +  llu-vlr  > 


(5.8) 


where  Cj.  is  a  constant  and 


inf  J  S  Hu.-uJI  X  (du  du9) 

^(CrC2)  f  "P 

0 


T 

Here,  *<(Cj,C2)  =  the  set  of  all  probability  measures  X  on  C^ 
prescribed  marginals  f.  and  Besides,  X  =  Xir  *  . 

1  «  v  l 


Ks(frf2)  = 


if  *  h 

if  q  =  C2' 

T 

,  x  C.,  with  the 


Proof A  complete  and  detailed  proof  of  this  result  will  appear  in  Baldwin,  et 
al.  [1].  Here  we  will  briefly  outline  the  basic  ideas  of  their  proof  with 
modifications  to  suit  our  needs. 

Let  A  =  (Sf(Y^):  X  €  ir(C^,)}.  Then  A  is  a  tight  subset  of  ir(c£.)  by 
Theorem  (2.4.3).  Define  the  map  '/'•tt(C^,)  — >  A  by  '/'(p)  =  Sf(Y^).  Again  by 
theorem  (2.4.3),  ^  is  sequentially  continuous. 
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T 

Note  that  A  is  a  subset  of  irfC^  ).  By  using  Theorem  (2.4.2)  for  k=l  and 

-P 

a  result  of  Mitoma  [14],  we  know  that  there  exists  an  index  p  £  p  such  that  A 

T 

is  tight  as  a  subset  of  ir(C^  _) .  Let  c£(coA)  denote  the  closure  of  the  convex 

-P 

T  T 

hull  of  A  in  t r(C^  _  ).  Then  c€(coA)  is  tight  in  ir(C^  _  )  and  is  therefore 

-p  ~P 

T  T  T 

compact  in  ir(C^  _  ).  The  canonical  inclusion  j:'w(C^)  _  )  C  *(c±-)  is 

-p  -P 

continuous  as  can  be  seen  from  the  proof  of  part  (c)  of  Theorem  6.3. 

T 

Therefore,  c£(coA)  can  be  viewed  as  a  subset  of  ir{ C^.)  and  is  a  compact  and 

T 

tight  subset  of  irfC^,). 

T 

It  can  be  shown  that  the  topology  of  weak  convergence  in  ir(C. , )  when 

T 

relativized  to  a  compact  tight  subset  of  Tr(C^,)  is  metrizable  so  that  c£(coA) 

is  a  Polish  space  under  this  topology. 

Let  p-  c£(coA)  — »  c£(coA)  be  the  restriction  of  p  to  c£(coA). 

c^(coA)  is  metrizable  and  so  sequential  continuity  of  p  is  equivalent  to 

continuity  of  p.  An  application  of  the  Schauder-Tychonof f  fixed  point  theorem 

T 

(see  [3])  gives  us  the  existence  of  the  McKean-Vlasov  equation.  If  €  7r(C^.) 
A 

such  that  A  =Sf(Y  ),  then,  for  this  choice  of  measure  A  ,  the  coefficients  A 
o  o 

and  B  satisfy  the  conditions  of  existence  and  uniqueness  of  solutions  as  listed 

in  Section  2.  Therefore  €  cj  ,  since  A  C  ir(c£  ). 

~P  T 

12  T 

For  part  (ii)  of  the  theorem,  let  A  and  A  be  two  measures  in  irfC^,)  such 
1  2 

that  \l=X(YX  )  and  A2  =  ^(Y*  ).  Then  Theorem  2.4.2  implies  that 

12  2 
E  sup  IIY*  II2  <  00  since  E  sup  IIY*  II2  <  “.  Likewise,  E  sup  IIY*  II2  <  ®. 
0£t£T  ~P  0£sfT  t  P  0£t<T  ~P 

A1  A2  2 

Therefore  if  Y  =  Y^  -Y  ,  E  sup  HY  II  <  00 .  Applying  the  I  to  lemma  to 
C  1  1  OisiT  P 

2  12 

IIYt.ll_p*exp(-2Gj,t)  where  Gj.  =  Gp(A  ,A  )  is  the  constant  that  appears  in  (5.8), 
and  then  using  the  condition  in  part  (ii)  of  this  theorem,  we  get  that 
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-T1 

EIIY^II_pexp(-2Crt)  <  E/J  KJX1 . X2) IIYsll_pGre  ‘  ds. 

Since  K  (X*.X2)  i  EIIY  II  .  the  above  inequality  yields 
s  s  -“p 

2  t  ”^ts  2 

EIIYpl_pexp(-2Cr)  Z  GjJq  e  J  EIIYsirpds. 

Gronwall’s  lemma  now  yields  EIIY2II  =  0  for  all  t  €  [O.T].  Since  Y  is  sample 

t  -p  t 

A1  X2  2 

continuous,  sup  IIYV  -Y  II  =0  a.s. 

CKtST  '  '  -» 

Remark  5.1:  Since  the  conditions  of  Theorem  5.1  hold  for  all  sufficiently 

large  m  with  p  being  the  index  determined  by  each  such  m,  the  conclusion  of  the 

theorem  holds  in  particular  when  m  is  replaced  by  a  larger  index.  Therefore, 

the  measure  Xq  obtained  in  Theorem  5.1  is  the  unique  solution  of  the 

McKean-Vlasov  equation  defined  by  equation  (5.3)  even  among  the  measures  in  the 

T 

larger  space  viz.  tt(C.  )  for  any  k  £  p.  This  fact  will  be  used  in  Section  6 

-k 

for  the  particular  choice  of  k=q  where  q  is  the  index  that  appears  in  Theorem 
4.1. 

§6.  PROPAGATION  OF  CHAOS 

T  ®n 

Let  denote  the  unique  probability  measure  on  (C^  )  that  solves  the 

B  -q 

martingale  problem  posed  by  the  system  of  equations  (3.2.1)  subject  to  the 

conditions  listed  in  subsection  3.3,  and  conditions  (SJC) ’  and  (MCir) .  Besides, 

T 

we  assume  that  b  and  c  satisfy  conditions  (5.6)  and  (5.7).  tj  €  Tr(ir(C^  ))  is 

-q 

given  by 

vB>  -  vf :  £  ,v " €  B)  v  B  €  ” 

i=l  yi  -q 

where  yn  denotes  a  generic  point  in  (C^  J®11  so  that  yn  =  (y^.y^. • • • .y^)  where 

-q 


in  proving  the  following  theorem. 


Theorem  6.1:  Under  the  conditions  specified  in  Theorem  5.1,  let  X  €  ir(C ^  )  be 

-q 

the  unique  probability  measure  that  solves  the  McKean-Vlasov  equation  (5.3). 

Then,  the  subsequence  (77  }  obtained  by  Theorem  4.1  is  such  that  rj  =>  6, 

\  \  o 

where  8  refers  to  the  Dirac  8  measure  provided  that  there  exists  a  >  0  such 
that  Elix”(0)ll^a  £  C,  where  C  is  independent  of  n. 

2  ~ 

Proof:  Let  f  €  (see  subsection  2.4)  so  that  f(u)  =  f(u[<p])  for  some 

<*>  €  <i>  and  f  €  C^(IR). 

Lij(f.yn.s)  =  f  '(yj(s)M){a(s.y”(s))[>]  +  b(s.yj(s),y”(s))[>]} 

+  £"(yj  (*)M) 

*  <J(j  1  (s)>  +  c(s.y"(s).y£(s)))*[.#>]. 

"  k2i^S,yj^8^  +  C^S’yj  ^S^’Yk 

where  yn  =  (y?.  .  .  .  .y11)  e  4>‘ x.  . . x4>' .  Let 
/v  a  n 


L(f.y(s).s.X)  =  f  ff’(y(s)[^]) 

T  L 


"<J> 


-q 


(a(s.y(s})(>]  +  b(s,y(s),z(s))[>]}  +  ^"(y(s)[>]) 


Q((X  (a(s.y(s))  +  c(s,y(s).z1(s))dz(z1))  J>]), 


(/  (<7(s,y(s))  +  c(s,y(s),z2(s))dX(z2))*[>])  dX(z) 


By  the  conditions  listed  in  subsection  3.3,  the  existence  of  a  unique  solution 

to  the  martingale  problem  posed  by  (3.2.1)  is  guaranteed.  We  have  called  such 

a  solution  as  P  .  Therefore 
n 

f(yj(t)[«Pl)  -  f(yj(r)0])  -  Uj  jCf.^.sids 

is  a  P^-martingale  with  0£r<t£T.  By  the  conditions  listed  in  the  statement  of 
Theorem  5.1,  a  unique  solution  of  the  McKean-Vlasov  equation  posed  by  (5.3) 
exists  and  is  denoted  by  Xq.  Therefore  the  following  is  a  Xo~mart ingale: 

f(y(t)0])  -  f((y(r)|>])  -  X\(f,y(s),s.Xo)ds 

T 

where  O^r^t^T.  Consider  the  function  F‘-  ir(Cp  )  — ►  IR  defined  by 

-q 

F(X)  =  s  (f(y(t)M)  -  f  (y(r)[f])  -  A(f  ,y(s)  ,s,X)ds}  (6.1) 

'T'  A 


g^yfr^).  ..g  (y(r  ))dX(y) 


1'  '  P  - '  7 

where  O^r, £r0.  .  . £r  £t<T  and  g, . g  are  bounded  functions  from  $  — *  IR. 

12  p  1  p  -q 

Hence  F(X  )=0.  Now  we  will  show  that  f  F^(X)q(dX)  =  0  by  direct 
o  T 

-q 

evaluation.  From  this,  it  follows  that  the  support  of  q  is  contained  in  the 
set  of  solutions  to  the  L-martingale  problem.  Corresponding  to  each  solution 
of  the  L-martingale  problem,  we  can  construct  a  weak  solution  of  the 
McKean-Vlasov  equation,  by  the  method  employed  by  Kallianpur  et  al.  [8].  From 
the  previous  section  we  know  that  the  McKean-Vlasov  equation  has  a  unique 
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strong  solution  namely  Y  .  Therefore  the  set  of  solutions  to  the  L-martingale 

problem  is  the  singleton  set  {Xq} .  From  the  fact  that  Xq  is  the  unique 

solution  to  the  martingale  problem  it  will  then  follow  that  tj=6^  . 

o 


Claim: 


:  lim  S  F^XJr)  (dX)  =  J  F2(X)t)(dX) 


*(<£  ) 

-q 


Proof  of  Claim.  Let  X  denote  Xt r  .  If  X  is  in  the  support  of  tj  recall  that 

s  s  nk 

T 

X  has  support  in  <i>_  for  each  0<s<T.  For  each  u,v  G  <P‘  ,  <p€<t>,  X  G  7r(Cp  )  and 

5  -q 


0<s^T,  let  for  R  >  0, 


a  (s,u)[<*>]  =  (-RVa(s,u)[if])AR 


b  (s,u)[<#>]  =  (-RVb(s,u.v))[<#>])AR 


Qh(s,u.X  )<*•*>  *  EAS>(s.u,XJ(f^) 
s  s 

~  R  R  R 

where  h(s,u,Xg)  =  ct(s,u)  +  c(s,u,Xs).  Replace  a.b  and  Q  by  a  ,b  and  Q 

respectively  in  the  definition  of  F,  and  call  the  resulting  function  as  F^. 


lim  f  F?(X)ir  (dX)  =  f  F?(X)T](dX) 
k-*»  _,„T  *  "k  ,  K 


^  w(cl  } 


) 


since  F^  G  C^irfC^  )).  The  claim  will  be  proved  if  we  show  that 


S  (F2(X)  -  F?(X))n  (dX))  and  S  (F^X)  -  F*(X))T}(dX)  can  both  be  made 

x  K  Hi,  T  K 

»(cj  )  »(cj  ) 

-q  -q 


arbitrarily  small  when  R  is  sufficiently  large. 


Using  Fubini's  theorem  and  Jensen's  inequality. 


S  (F2(A)-f2(X))ti  (dX)$3(t-r)(  max  llg  ll2)P 
-q 
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X  JX  C  X  (f(y[*]a(s.y)[*])\(dy)  (6.2) 

r  X  {y:|a(s.y)[^]|>R} 


{y:  |b(s.y,Xs)[«f]i>R} 


(f(yW)b(s.y.X,)W)\(dy) 


*  ,  „  !  ,  4  («"(yM)«h(s,y,x  Vdy»ds 

{y;Qh(S.y.Xs)(,''f)>R> 

It  remains  to  show  that  each  of  the  three  terms  on  the  right  side  of  (6.2)  can 
be  made  arbitrarily  small,  uniformly  in  k,  when  R  is  large.  Since  the  method 
for  each  term  is  essentially  the  same,  we  shall  consider  only  the  third  term. 

V.y.X  *  M'm  |Qh(s.y.X  )U.-»  *  M&3e( lMllyll?n|+K2(n,))+3C2(m)] 


by  using  (SLG)  and  equation  (5.7).  Thus 


<y:  Si(s,y,X  )(»•»)  >  R>  e  (y;  »y|Cm  >  W 


where  k  is  a  suitable  positive  constant  and  R  is  sufficiently  large. 
Therefore,  the  third  term  on  the  right  side  of  (6.2)  is 


*  K:  sl  S  .  .{  „2  <ullyll-«)Vdy)ds 

\  {y:l!yll_m>R/k> 


by  using  (SLG)  again  with  as  a  suitable  constant  independent  of  t^. 
Continuing: 

*  K1  K  t  £  ^  ,  Jds 

K  J=i  { llx j  ( s )  ll_m  >  R/k) 

S  K.24+a  ~  ^  (E(l+llxI)k(s)ll4+a)}4/^4+a^{P(llxI!k(s)ll2  >  R/k)}a/^4+a^ds 

*  1  r  n,  1  v  j  v  /  -m  11  1  v  j  ’  -m 


by  Holder’s  inequality. 
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i  k3/Sa/‘4*a)  —  0 


as  R  — *  00 , 


wherein  the  last  inequality  uses  Chebyshev’s  inequality  and  the  moment  bound 
given  by  (3.3.1).  Besides,  is  a  constant  independent  of  n^.  The  fact  that 
X  (F^A)  -  F?(A))77(dA)  can  be  made  small  for  large  R.  follows  along  the 

rp  K 

*(<£  ) 

-q 

same  lines  as  above  and  its  proof  is  hence  omitted. 


Continuation  of  the  proof  of  the  theorem:  Using  the  above  claim, 

sr  F(A)2Tj(dA)  =  lim  X  «n.  [A  .^(f  {y\  t)[>]) 
ir( Ci  )  k-*»  ir(C^  )  K  "k  J 


-q 


-  7(y^(r))W  -  /'  ij-  ,j(f.A.)J.}e1(y^(r1))...Bp(y^k(rp))]2dPllk(ynk) 
since  rj  (B)  =  P  (y^:  - —  2^  6  €  B)  so  that 

V  v~  "k  i=i  "k 

yi 

x  ^  -  f(y?(r)[^])  -  xj  i-  (f.y^.sjds) 

t(cT  J  j  rnkl-'1J  - 

-q 

•  g1(y^(r1))...gp(y^(rp)))%  (i-  £  6  ) 

yi 

%  ((f(y?(OM)  -  fCy^rJC-p])  -  sl~  (f.y^.sjds) 

nf  J  J  r  \  i=l 


=  X 


TT(CJ  )®”k  "k 
-q 


*i(y^(rI))...gp(y^(rp)))2dpnk(ynk) 


(6.2) 


Then,  the  right  side  of  (6.2)  is  equal  to 


S  /!  4  »g.(y"k(r  )). 

y(CJ  )%r\  ' 

-q 


.gp(y^k(rp)))2((f,(y^k(s)M))2 
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*  Q((Jr  ^{^(s.y'Xs))  +  c(s,yI!k(s).y^k(s))})i<[^], 
i=l  J  J  1 

(~  ^{^(s.y^Cs))  +  cCs.y^Cs)  ,y^k(s))})M[<p])dsdP  (y1^)  =  o(^r) 

\  1=1  J  J  1  nT 


by  using  the  condition  (SLG)  and  the  bound  given  by  (3.3.1).  Thus 


J*  F(X)2T)(dX)  =  lim  (o(i-)  +  X 

oT  x  kn»  "k  \  ,„T 


»(S  ) 


(cj  i®"* 


(f(yx  (tJW) 


-  f(y^(r)W)  -  si  ~  ^  L*  (f.y^.sjds) 

1  rnki=l  1,1  w 

•  (f(y^(t)W)  -  f(y^(r)M)  -  X*  (f.y^.sjds) 

Z  z  r  \  i=i  Xt*  ~ 

gl(y?{rl))*  •  •«p(ylk(rp))gl(y2k(rl))*  •  •gp(y2k(rp)) 

dP  (y"14)}  =  0 
"k  ~ 

since  the  independence  of  the  Wiener  processes  W*  and  W2  implies  that 

<M*(t),  M^(t)>  =  0 

where 

M*(t)  =  f(y^(t)M)  -  f(y^(r)[v>])  -  X*  ^  2^  L* ^(f .y^.sjds 


*J(t) 


which  is  a  P  -martingale. 

"k 

Thus  X  F(X)277(dA)  =  0  for  all  F  defined  by  (6.1)  with  f  €  2J2(4>'), 

T  D 

*(S  ) 

-q 

p  €  IN  and  gj . gp  continuous  and  bounded  mapping  <£'  — *  1R,  and 

T 

0£rj£...£r  it$T.  Since  Xq  is  the  unique  member  of  tt( C ^  )  such  that  F(Xq)  =  0, 

p  -q 
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we  get  that  q  =  5*  . 


Remark  6.1.  By  Theorem  6.1,  the  possibly  random  measure  p  such  that 

l 

"k 


D  T 

— ►  H  €  has  been  shown  to  be  a  non-random  measure  and  in  fact 


Our  result  on  propagation  of  chaos  is  presented  in  the  next  theorem. 


Theorem  6.2.  Under  the  conditions  (SA) ,  (SIC),  (SCC),  (SLG),  (SJC),  (SMC), 

(SJC) '  and  (MCir)  and  with  the  coefficients  b  and  c  satisfying  (5.6)  and  (5.7), 

and  sup  Elix"(0)ll^+a  i  C  for  some  a  >  0  we  get 
n  ' 

(i)  rjn  =>  in  w(w( )). 

0  -q 

(ii)  If  T]t  =  iffp*")  where  pt(w,Xn)  =  —  2  I 

n  n  n  ~  n  i=l  X"(t.w) 


and  (X?(*),  l£i<n)  solves  the  system  (3.2.1),  then  q*  =>  6,  ir  *  for  0<t£T. 
i  n  Aq  t 

That  is  to  say  p*  — »  X^tt^  *n  distribution  hence,  in  probability  as  well. 
Proof.:  We  have  shown  in  Theorem  6.1  that  v  =>  6,  .  In  fact,  for  anv 

\  xo 

convergent  subsequence  (p  }  of  the  sequence  of  empirical  measures  {p  },  we 

V  n 

get  from  Theorem  6.1  that  q  =>  6^  .  Therefore,  the  whole  sequence  q  weakly 

n j 1  xo  n 

converges  to  6,  . 

\> 

To  prove  (ii),  note  that  for  all  real-valued  continuous,  bounded  functions 
f  on  tt(cJ.). 


f  f(X)qn(dX) 


-q 


*(c l  ) 

-q 


f(X)5  (dX). 

Ao 


(6.2) 
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In  particular,  if  f(A)  =  S  g(y)dA(y)  where  for  all  y  €  ,  g(y)  =  g(y  ) 

(cj  )  ~q 

-q 

for  t  fixed  in  [O.T],  and  g  a  continuous  bounded  function  from  <t>_  to  IR,  then 

T 

f(A)  is  indeed  a  real-valued,  continuous,  bounded  function  on  C^,). 

Therefore,  (6.2)  for  this  choice  of  f  implies  that  r?1"  =>  5,  it  *.  Thus 

n  XQ  t 

p^  — »  AqV^  aS  n  — *  °°  *n  probability  since  the  limit  is  non-random. 


n 


Theorem  6.3:  (a)  For  each  T>0,  let  the  conditions  (SA) ,  (SIC),  (SCC) ,  (SLG) , 
(SJC)  and  (SMC)  hold.  Then  the  system  of  SDE's  (3.2.1)  admits  a  weak  solution 
that  is  pathwise  unique.  That  is,  (3.2.1)  has  a  unique  strong  solution  in 

(C„  )*”■ 

(b)  Assume  the  additional  conditions  (SJC)'  and  (MCtt)  for  each  T>0.  If  b  and  c 
are  as  specified  by  (5.6)  and  (5.7),  the  McKean-Vlasov  equation  posed  by  (5.3) 
has  a  unique  strong  solution  in  C^,. 

1  n 

v“  T  so  that  its 


(c)  In  view  of  (a)  and  (b)  above,  define  p  (*.w)  =  —  21 

n  n  i=l  Xj(* ,w) 

law  t}^  €  ir(v( C^,))  n  £  1.  Let  Xq  be  the  probability  measure  on  C^.  that  solves 

the  McKean-Vlasov  equation  posed  by  (5.3).  If  sup  Elix!j1(0)ll_*a  £  C  for  some 

n 

a  >  0,  then 


Vn  =>  ln  ir(ir(C4,-))' 


Proof :  Part  (a)  follows  by  reading  off  the  corresponding  result  in  Kallianpur, 
Mitoma  and  Wolpert  [8]. 

(b)  Let  the  conditions  of  Theorem  6.2  hold  for  each  fixed  T>0.  Then,  the 
results  in  Section  5  and  6  hold  in  the  interval  [0,®).  To  see  this,  suppose 
\) 

Xq  =  S£(Y  )  solves  the  McKean-Vlasov  equation  in  the  interval  [0,Tq],  and 
-X1 

Aj  =  !£(Y  )  solves  the  McKean-Vlasov  equation  in  the  interval  [0,T],  where 
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Tj>Tq,  then,  by  the  uniqueness  of  solutions  to  the  McKean-Vlasov  equation,  we 
get  that  the  projection  of  on  the  interval  [O.Tq]  must  coincide  with  Xq. 
Thus  Xj  is  an  extension  of  Xq  in  the  above  sense.  Such  an  argument  shows  the 
existence  and  uniqueness  of  solutions  to  the  McKean-Vlasov  equation  (5.3)  in 
the  interval  [0,T]  for  any  T>0. 

Choose  Tn  =  n,  and  the  corresponding  measures  X^e^C^,).  Solve  the 
McKean-Vlasov  equation  (5.3).  Then,  the  projective  limit  of  {X^}  is  a  measure 
Xq  €  7r(Cp,)  that  solves  (5.3)  for  all  t  >  0. 

To  prove  part  (c),  we  make  the  following  observations:  For  any  positive 

T  T  T 

index  k,  i:  C  C^,  is  continuous.  In  fact,  the  topology  on  C^,  as  given 

in  Section  2  is  equivalent  to  the  weakest  topology  with  respect  to  which  the 
above  canonical  inclusions  are  continuous.  Therefore,  we  have 


Claim:  If  f  =>  r  in  w(tt(cJ  )).  then  T  =>  T  in  ir(w(cj,)). 

-k 

T  T 

Proof  of  the  claim'-  First,  note  that  the  inclusion  j:  w(C^  )  C  tt(Ca.)  defined 


-k 


$>' 

-1  T 

by  j(X)  =  Xi  is  continuous.  To  see  this,  let  (X  }  be  a  net  in  ir(C.  )  such 

a  -k 

T 

that  X  =>  X  in  w(C,  ).  Therefore 


-k 


J  fd*«  -  Sr  fdAVf  CChtcJ  ). 

^T  JT  -k 


-k 


-k 


T  T 

Let  g  €  (^(Cp,).  The  composition  g*i  is  then  in  (^(C^  )•  Also 

-k 

X  gdX  i  *  =  X  g*idX  for  all  a 
_  a  T  a 

C*  ci 

<J>  <P  , 

-k 


and 
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so  that 


f  gdX  =  J  g*idX 

T  T 

c1  c1 

tf>  <J> 


J  gdXa  —  s  gdX  for  all  g  €  <yC^.). 


Therefore  j  is  continuous. 

T  T 

Now  let  k  be  the  inclusion  from  t(tt(C.  ))  C  tr(ir(C.  ,  )) .  Continuity  of  k 

-k 

can  be  proved  by  following  step  by  step  the  proof  of  continuity  of  j.  The 
claim  is  thus  shown. 

Now  part  (c)  is  shown  by  observing  that  Theorem  6.2  part  (i)  implies  that 


T]  =>  6.  in  ir(ir(C*  ) 
n  A0  % 

and  hence 


by  the  claim  shown  above. 


^n  =>  \  in 


Vn=>  6\q  in  T'(T(c<f) 


since  the  inclusion  ) )  C  ) )  is  also  continuous.  Note  that  r?n 

—  T  — 

and  6,  are  the  projections  of  tj  and  5r-  on  7r(Tr(C.,)).  Thus  tj  =>  6r-  in 
X0  n  A0  *  n  A0 

Remark  6.3:  The  unique  strong  solutions  mentioned  in  parts  (a)  and  (b)  of  the 
above  theorem  are  in  general  <P‘ -valued  processes  and  cannot  be  guaranteed  to 
lie  in  a  single  Hilbert  space  .  This  is  so  since  the  indices  m  and  p  vary 

«J 

with  T  in  the  conditions. 


§7.  APPLICATION  TO  INTERACTING  SYSTEMS  OF  NEURONS. 


The  random  behavior  of  the  voltage  potential  of  a  spatially  distributed 
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neuron  has  attracted  considerable  attention  in  neurophysiology  and  can  be 
modeled  in  the  following  set-up: 

Let  H  be  a  separable  Hilbert  space  and  T  be  a  strongly  continuous 

contraction  semigroup  on  H  with  a  densely  defined,  closed,  negative-definite 

2 

generator  A.  In  practice,  H  is  usually  taken  to  be  L  (3t,p)  where  3C  is  the 

membrane  of  a  neuron  and  p  is  a  suitable  measure  on  3t.  If  there  exists  r^>0  so 

-rl  oo 

that  (I-A)  is  Hilbert-Schmidt,  then  there  exists  a  for  H  such 

oo  “2ri 

that  -  A<#> .  =  A  <p  ,  j=l,2,...  with  2  (1+A.)  <  ®°.  Let 

J  J  J  J 


❖  =  {<p€H:  2  (1+A  )^r(<p,<p  .)u  <  “  for  any  r>0}.  Define  on  $  a  family  of 

j=l  J  J 


increasing  Hilbertian  norms  ll»ll  with  ll<pll 


?  =  2  (l+A.)2r(-p,*.)2  and  let  ❖ 

j=l  J  J 


denote  the  completion  of  ❖  w.r.t.  If •  II  .  Since  ❖  C  ❖  is  Hilbert-Schmidt, 

r  r+r j  r 

it  is  easy  to  see  that  Ms  a  nuclear  space.  The  semigroup  {Tt}t^o  can  be 
written  as  follows.  For  any  <#>  €  ❖ 


00 


V  * 


e  ❖. 


The  voltage  potential  is  identified  as  the  solution  of  the  ❖’-valued  SDE. 


dX  =  A’Xtdt  +  dW 

where  A’  is  the  adjoint  operator  of  A  and  W  is  a  ❖’-Brownian  motion  with  a 
certain  covariance  function  E(Wt,<#>)(Wj  ,❖)  =  (t  A  s)Q(« p.'p). 

More  generally,  suppose  At  generates  a  strongly  continuous  contraction 
evolution  operators  T(s,t),  s£t  on  ❖.  Assume  the  following  conditions  on  At= 
For  any  T  and  large  enough  m,  there  exists  a  p>m  such  that  |A^.u|_p<  K|u|_m  for 
all  t£T,  and  u€$’ ,  i.e.,  as  continuous  linear  operators  from  ❖’  to  ❖’ .  (A  }  /T 

m  in  p  l  l  s  i 

are  uniformly  bounded. 
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Then  the  following  <J>’-valued  SDE  modeling  the  voltage  potential  of  a 
neuron  has  a  unique  solution: 

dXt  =  AtXtdt  +  dWt 

=  f  . 

0  4o 

Moreover,  the  solution  can  be  explicitly  written  as: 


xt  =  r(o..)f0  +  »t*j'Asr,twsds. 


Here,  T'  denotes  the  adjoint  operator  of  T  .  Now,  consider  the  system  of 

S  |  t  S  i  C 

n-interacting  neurons  whose  voltage  potentials  are  governed  by  the  following 

SDE: 


dx[n)(t)  =  (Atx[n)(t) 


Xj(t))dt 


dW^t),  i=l  ,2, 


,n  (7.1) 


Xn(0)  =  f  €  <P'  where  b  :  <P'x<P'  — ►  <£'  represents  the  interaction  between  neurons 
and  {W^t)}^^  are  independent  copies  of  a  ^'-valued  Brownian  motion. 

We  require  that  the  interaction  b^.:  satisfy  the  conditions 

(SOC),  (SLG) ,  (SJC),  (SMC).  (SJC) '  and  (MCrr)  as  in  Theorem  6.2. 

The  existence  and  uniqueness  theorem  in  Section  3  thus  guarantees  that 
system  (7.1)  has  a  unique  solution.  The  propagation  of  chaos  in  Section  6 

1  n  y 

asserts  that  the  empirical  distribution  —  2  5  ,  €  ir(C.  )  converges  in 

n  i=l  Xjnj  -q 

T 

probability  to  a  deterministic  probability  measure  Xq  €  ir(C^  )  which  is  the 

-q 

law  of  the  solution  of  the  McKean-Vlasov  equation  corresponding  to  (7.1): 


dXt  =  (AtXt  +  bt[Xt,xJ])dt  +  dWt, 


b[X,xh  =  f  b(x.y)dxhy) 

<P‘ 


where 
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and  Aq  is  the  law  of  X{. 

Thus  the  asymptotic  behavior  of  a  large  system  of  neurons  through 
mean-interactions  becomes  asymptotically  independent  with  the  distribution 
governed  by  the  McKean-Vlasov  equation  (7.2). 

8.  CHAOTIC  SYSTEMS 

8. 1  Exchangeable  systems. 

Till  now,  the  initial  random  variables  Xj(0) .  l£j<n  have  been  assumed  to 

be  i.i.d.  random  variables.  We  now  relax  this  condition  and  assume  that  X^(0), 

l£j<[n  are  exchangeable  random  variables  for  each  n£l.  That  is,  the  law  of 
n  n  &n 

Xj(0)  l<j<n,  denoted  by  Pq  €  v(9')  is  a  symmetric  probability  measure  on 

We  call  the  symmetric  measures  p^  p^-chaotic  if  the  following  condition 
holds:  For  every  integer  k£l  and  fj . fk  €  Cb(<f), 

**"  . fk(uk)d4o(“>  *  {  WW  <81> 

n-*°  (9  )  i=l  9 

where  u  =(u^ . u^)  €  ( 9 ')  ,  and  Pq  is  a  probability  measure  on  9' .  We 

assume  that  the  measures  p^  are  p^-chaotic . 

In  the  context  of  the  neuorphysio logical  model  described  in  section  7,  the 

assumption  of  exchangeability  of  the  law  of  (X^(0) . X^(0))  for  each  n^l  is 

equivalent  to  saying  that  the  particular  order  in  which  the  neuronal  membranes 
are  taken,  is  immaterial.  This  is  so  since  the  random  variables 

Xj(t) . X^j(t)  eac^1  *20  n£l  turn  out  to  be  exchangeable  random 

variables.  The  pQ-chaoticity  assumption  is  needed  in  showing  the  propagation 
of  chaos  result.  To  see  this,  consider  the  simplest  case  where  the  drift  and 
diffusion  coefficients  are  identically  zero  so  that  (8.1)  itself  becomes  the 
propagation  of  chaos  statement. 
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The  results  of  the  previous  sections  hold  for  the  exchangeable  model  as 
well  if  we  assume  Pg-chaoticity  and  that  EIIX^(0)ll^*d  £  C  where  C  is  a  constant 
independent  of  n,  and  m  is  the  index  that  appears  in  the  conditions  listed  in 
subsection  3.3. 

8.2  Finite-dimensional  systems 

d  2^2  d 

By  setting  $  =  1R  with  llxll  =  2  xf  for  each  n^l,  <t>  =  Dr  =  D  <t>  is  seen 

"  1=1  1  nil  n 

to  be  a  nuclear  space  with  its  strong  dual  being  isomorphic  to  IRd.  In  this 
case,  all  the  norms  II  *11^.,  -40  <  k  <  ®  are  one  and  the  same,  namely,  the 
Euclidean  norm  on  IRd  denoted  by  11*11.  Therefore,  the  indices  m,p,q  etc.  in  our 
conditions  can  and  will  be  taken  to  be  1.  The  canonical  maps  j  will  not 
appear  in  the  conditions  in  this  case.  Besides,  expressions  such  as 
|Q  /t.  ,,-v  l_m  __  will  simply  read  as  trace(<7CT*(t,u)) .  Also,  the  condition  (SA) 

\J  (  t  |  Uy  — III  i  ill 

J 

is  trivially  seen  to  hold  for  the  choice  of  $  =  R.  The  propagation  of  chaos 
result  for  the  finite-dimensional  exchangeable  system  is  given  in  the  next 
theorem: 

Theorem  8.2.1.  For  each  T  >  0,  let  the  conditions  (SIC),  (SCC) ,  (SLG) ,  (SJC) 
and  (SMC)  hold.  Let  Xj(0)  l£j£n  be  exchangeable  random  variables  and  let 
p”  =  law  of  (X^(0), — ,X^(0))  be  |iQ-chaotic.  Then, 

(9l  - 

a)  The  system  of  SDE’s  (3.2.1)  has  a  unique  strong  solution  in  (C  ,) 

nr 

b)  In  addition,  assume  (SJC)'  and  (MCjr)  for  each  T  >  0.  If  the  coefficients  b 

and  c  are  as  specified  by  (5.6)  and  (5.7),  then  the  McKean-Vlasov  equation 

(5.3)  has  a  unique  strong  solution  in  C  ,. 

Rd 

c)  Assume  the  conditions  in  part  (b).  Further,  assume  that  there  exists  6>0 
such  that  EIIX^(0)ll^+d  £  C  where  C  is  independent  of  n.  Then,  in  the  notation 
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of  section  6 

17  =>  5  in  7r(rr(C  ,)) . 

"  X  IRd 

o 

The  above  result  enables  us  to  compare  our  results  with  those  of  Sznitman  [16]. 
The  conditions  made  by  Sznitman  are  the  following: 

(i)  The  initial  random  variables  X^(0) . X^(0)  are  IR^-valued  exchangeable 

random  variables  and  are  bounded. 

(ii)  p”  =  law  of  (X^(0) . X^(0))  on  (IR^)®11  are  ^-chaotic,  where  is  a 

probability  measure  on  IR^ 

(iii)  The  drift  and  diffusion  coefficients  are  uniformly  bounded  and  satisfy 
uniform  Lipschitz  conditions  in  the  space  and  time  variables. 

(iv)  The  covariance  form  Q  is  the  identity  matrix. 

In  the  next  paragraph  the  conditions  (i)  through  (iv)  are  compared  with  those 
that  appear  in  Theorem  8.2.1. 

First  (SIC)  and  the  moment  condition  introduced  in  part  (c)  of  Theorem 
8.2.1  are  satisfied  since  (i)  says  that  the  initial  variables  are  bounded. 

(SOC)  is  verified  as  follows: 

For  u.v  €  IR^,  O^t^T,  and  h(t,u,v)  =  a(t,u)  +  c(t,u,v), 

|2a(t,u)*u  +  2b(t,u,v)*u  +  tr(hh*(t,u,v)) | 

£  2llalljlull  +  Sllbll^llull  +  trace(hh*(t,u,v)) 

<,  0(  1+llull2) 

by  using  the  uniform  boundedness  of  the  coefficients  and  condition  (iv). 

The  verifications  of  (SLG)  and  (SJC)  given  the  Conditions  (i)  to  (iv)  are 
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simple  and  hence  left  to  the  reader.  (SMC)  can  be  verified  by  using  Lipschitz 
continuity  of  the  coefficients.  (SJC)  and  (MCir)  are  obtained  by  continuity  and 
boundedness  of  the  coefficients.  Thus,  our  set  of  conditions  for  the 
propagation  of  chaos  is  weaker  than  that  imposed  in  the  finite-dimensional 
set-up  by  Sznitman  [16].  The  finite-dimensional  result  of  Leonard  [11]  is 
close  in  spirit  to  Theorem  8.2.1  and  hence,  a  comparison  of  the  two  is  left  to 
the  reader. 
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